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This paper reexamines the physical roles of trapped and passing electrons in Bernstein-Greene- 
Kruskal (BGK) electron solitary waves, also called the BGK electron holes (EH). It is shown that 
the shielding of the positive core is achieved by electrons trapped and oscillating inside the potential 
energy trough, instead of the thermal screening from the ambient plasma as previously thought. 
The size of a BGK EH is therefore not restricted to be greater than the Debye length A_d. 



PACS numbers; 52.35.Sb, 52.35. Mw, 52.35.Fp 

In 1957, Bernstein, Greene and Kruskal solved the 
one-dimensional, time-independent Vlasov-Poisson equa- 
tions and obtained the general solutions for electrostatic 
nonlinear traveling waves, including solitary potential 
pulses. Their derivation emphasized the special role 
played by the particles trapped in the potential energy 
troughs. They demonstrated that one could construct 
waves of arbitrary shapes by assigning the distribution 
of trapped particles suitable for the desired wave form. 

In 1967, Roberts and Berk provided a quasi-particle 
picture for the electron phase space holes (EH) based 
on the results of a numerical experiment on two-stream 
instability. They followed the evolution of the electron 
phase space boundary between / = and / = 1 us- 
ing the time-dependent Vlasov-Poisson equations with 
the "water-bag" model in which the distribution function 
only takes discrete values, and 1. They interpreted the 
elliptical empty (/ — 0) region associated with a positive 
charge observed in the late stage of the nonlinear develop- 
ment as a BGK EH. In order to explain the coalescence of 
neighboring EHs, a negative effective mass was assigned 
to each EH to compensate for the Coulomb repulsion of 
two positive EHs. Thus, they suggested a quasi-particle 
picture that BGK EHs have positive charge and negative 
mass. This picture is in use even today to interpret re- 
sults in computer simulations of EH disruptions |^ , and 
to model the mutual interaction of electrostatic solitary 
waves in space plasma 

It was not until 1979 that BGK EHs were experimen- 
tally reahzed by the Ris0 laboratory experiments [§,1). 
By applying large amplitude potential pulses in a plasma- 
loaded wave guide, solitary potential pulses were excited, 
including EHs and Korteweg-de Vries (KdV) solitons . 
Investigations of the mutual interaction of EHs showed 
that two EHs close to each other would coalesce if they 
have almost equal velocity and they would pass through 
each other if their relative velocity was large §. The 
coalescence was interpreted in terms of the positive EH 
picture derived earlier |^,^ . 

The analytical work that followed the Ris0 experi- 
ments mainly focused on constructing solutions and ob- 
taining the corresponding width-amplitude relations to 



facilitate the comparison between BGK EHs and KdV 



solitons [§ 10|. In addressing the quasi-particle picture 
of EHs, Schamel j|] concluded that EHs were positively 
charged, screened by the ambient electrons over many 
Debye lengths (A_d), and had negative mass [||. This 
conclusion supports the positive EH picture previously 
obtained from the numerical experiments ||,|ll | and that 
the minimum size of EHs of several Ad is a consequence 
of thermal screening by the ambient electrons. 

Turikov [|| followed the BGK approach and con- 
structed the trapped electron distribution for a 
Maxwellian ambient electron distribution and various 
solitary potential profiles. He restricted his study to 
BGK EHs with phase space density zero in a hole center 
and the results showed that the potential width increases 
with increasing amplitude. This behavior is different 
from that of the KdV solitons whose width decreases with 
increasing amplitude. He also numerically simulated the 
temporal evolution of the EHs for different Mach num- 
bers to study the EH stability and found that EHs are 
quasi-stable for Mach numbers less than 2. One of the 
main conclusions that Turikov made was that EHs are 
purely kinetic nonlinear objects in which trapped elec- 
trons play an important role, but exactly what physical 
role trapped electrons played was not addressed. 

Space-borne experiments now show that electrostatic 
solitary waves (ESW) are ubiquitous in Earth's magne- 
tospheric boundaries, shock, geomagnetic tail and auro- 
ral ionosphere [p^-|l7| . While detailed properties of these 
solitary waves are continuously being studied, it has been 
shown that in a number of cases the ESWs have features 
that are consistent with BGK electron |l^ and ion |l^] 
mode solitary waves. A statistical study by FAST satel- 
lite observations |l^,|l^ in the auroral ionosphere has re- 
vealed that solitary pulses with a positive potential typ- 
ically have a Gaussian half width ranging from less than 
one Ad to several Ad with a mean of I.SOAd and a stan- 
dard deviation of 1.13Ad- However, the statistical anal- 
ysis strongly favored large amplitude pulses [|l5|,0 and 
smaller EHs if they existed were not sampled. This ques- 
tion of how small EHs can be is an important issue asso- 
ciated with how a coUisionless plasma supports nonlinear 
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waves and needs to be further investigated. 

To examine the physical roles played by trapped and 
passing electrons, we start first with simple physical ar- 
guments and then perform analytical calculations using 
the same formulation used by Turikov except that 
we relax the restriction of empty-centered EHs and ob- 
tain a more general width-amplitude relation. The num- 
ber density profiles of trapped and passing electrons are 
calculated to quantify the separate contributions from 
trapped and passing electrons to the charge density. Im- 
plications of our results to the positive EH picture and 
minimum size of BGK EHs are made. 

We first discuss heuristically the behavior of electrons 
in the vicinity of a potential pulse. Figures |l|(a)-|l|(c) 
show the general form of a positive solitary potential 
pulse {(j){x)), the corresponding bipolar electric field {E = 
—d(t)/dx) and the total charge density (p = —d^(j)/dx^). 
The charge density is positive at the core, negative at the 
boundary, and zero outside the solitary potential. Figure 
|l|(d) shows the potential energy trough with an electron 
passing by (open circle) and a trapped electron (solid cir- 
cle) at its turning point. Consider the phase space tra- 
jectories of electrons passing by the potential and those 
that are trapped in the potential shown in Figure |l|(e). 
The dashed line marked by electrons with zero total en- 
ergy is the boundary of the trapping region inside which 
electrons are trapped and outside which electrons are un- 
trapped. The total energy, w = — e(j), is a constant 
of motion. A passing electron {w > 0) moves with a 
constant velocity outside the potential and the speed in- 
creases when it encounters the potential pulse and then 
decreases back to its original value as it moves away. A 
trapped electron {w < 0) bounces back and forth between 
its two turning points in the potential. Since there is no 
source or sink for the particles, the density is inversely 
proportional to the velocity. We can thus deduce that the 
density of passing electrons is constant outside and be- 
comes smaller as (p increases. No excess negative charge 
results from the passing electrons. On the other hand, 
trapped electrons have density maxima at their turning 
points, and so must be responsible for the excess nega- 
tive charge. The charge density variation (Figure 0(c)) 
needed to be self-consistent with the solitary potential 
pulse is thus a net balance of the negative charge from 
trapped electrons and the positive charge due to the den- 
sity decrease of passing electrons since the ion density is 
assumed uniform. From this simple picture, one can see 
that in BGK solitary waves it is the trapped electrons 
traveling with the solitary potential that screen out the 
positive core. Our picture is different from the picture of 
a positive object in a plasma whose screening is achieved 
by the thermal motion of the plasma (Debye shielding). 

The entire trapping region that consists of the total 
electron density enhancement at the flanks and deple- 
tion at the core is a physical entity produced by the self- 
consistent interaction between the plasma particles and 



the solitary potential pulse. This defines the physical 
identity of one BGK EH. The total charge of the entire 
trapping region is zero, and therefore it follows that two 
separated BGK EHs do not interact and the concept of 
negative mass is not needed. 

We now use the approach formulated by BGK 0| to 
quantify the above arguments and further demonstrate 
that the results are independent of the strength of the 
nonlinearity. The time-independent, coupled Vlasov and 
Poisson equations with the assumption of a uniform neu- 
tralizing ion background take the following form: 
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(1) 



(2) 



where f is the electron distribution function and the quan- 
tities have been normalized with the units of the Debye 
length X]j, the ambient electron thermal energy Tg, and 
the electron thermal velocity vt = \/2Tf,/m. The total 
energy w = — (j) under this convention. Any / = f{w) 
is a solution to Eq. |l| as can be readily verified. Rec- 
ognizing this, Eq. |2| can be re-writtcn in the following 
form. 



dx"^ 



2^w + (f) 



dw 

+ (p 



-1, (3) 



where ftr(w) and fp(w) are the trapped and passing elec- 
tron phase space densities at energy w, respectively. The 
first integral on the RHS of Eq. ^ is the trapped elec- 
tron density, and the second integral the passing electron 
density. Prescribe the solitary potential as a Gaussian, 



(/)(V',5,x) = Vexp(-a;V2(5^), 



(4) 



and the passing electron distribution a Maxwellian where 
the density has been normalized to 1 outside the solitary 
potential. 



fp{w) = ^exp(-u;). 



(5) 



Following the BGK approach, we obtain the trapped elec- 
tron distribution. 



ftr{lp,S,w) 




(6) 



The first term in ftj. stems from d^4>/dx^ term in Eq. 
^ and has a single peak &t w = J^S72- This term is 
at w = 0^, goes negative at w — —■(/;, and will al- 
ways be single peaked even for other bell-shaped soli- 
tary potentials (for example, sech^{x/6) and sech^(a;/<5). 
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see Figure 4 of Q for the special case of empty-centered 
EHs). Although the peak location may vary, it will not 
be at the end points, and —-ip. The second term aris- 
ing from the integral of the passing electron distribu- 
tion decreases monotonically from w = 0~ to w = —ip. 
The end point behavior of the two terms implies that 
ftriw = 0^) > ftr{w — —ip). Combining the behav- 
ior of the two terms in ftr, it can be concluded that 
ftr{0 > w > — ?A) > ftr{w — This feature of ftr 

is essential in making a solitary pulse, and it manifests 
itself at the peak of the potential as two counterstream- 
ing beams. For the solution to be physical, ftr has to be 
nonnegative. ftr{'>jj,S,'w = —tp) > suffices this require- 
ment, yielding the width-amplitude (S-ip) relation, 



d > 



2V?(21n4- 1) 



1/2 



(7) 



The equal sign in Incq. ^ corresponds to the case of 
empty-centered EHs studied by Turikov We plot 
Ineq. |^ in Figure ^j. A point in the shaded region repre- 
sents an allowed EH with a given ip and S. The shaded 
region includes all of the allowed tp and 6 for the range 
of values shown. For a fixed (5, all < "00 are allowed, 
where ipQ is such that ftritpOy^T^ — — ?Ao) — 0; while 
for a fixed ip, all S > Sq are allowed, where Sq is such 
that /tr(V': ^OiW = —ip) = 0. The corresponding physical 
meaning is that empty-centered EHs give the largest am- 
plitude (ipo) for a fixed width, and the smallest width (Sq) 
for a fixed amplitude as comparing to EHs that are not 
empty-centered. Note that by distributing a finite num- 
ber of electrons at rest at the bottom of the potential 
energy trough, the plasma supports smaller amplitude 
or larger width pulses than the amplitude and width of 
empty-centered EHs. This inequality relation is dramat- 
ically different from that of KdV solitons whose width- 
amplitude relation is a one-to-one mapping. 

With fp and ftr, we can now calculate the passing and 
trapped electron densities separately and obtain 



np{x) = exp {(p) l-erf(\/^) 



ntrix) 



-0[l + 21n(0/7/;)] 

(52 



+ l-exp{(p) l-erf(V^) 



Another way to obtain ntr comes directly from solving 
Eq. H as the first term on the RHS is exactly ntr- Solv- 
ing Eq. ^ for ntr is simpler, but without the knowledge 
of ftr, one is not guaranteed whether the particular set 
of (-0, S) is physically allowed {ftr > 0). Note that even 
for < 1, n„ 1 - 2 



\/(p/n is different from the lead- 
ing terms of a Boltzmann distribution, e"^. The physical 
meaning is that under self-consistent interaction of elec- 
trons and the solitary potential, electrons in the vicinity 
of the potential are not in local thermal equilibrium, in 



contrast to the starting point of local thermal equilibrium 
in obtaining the thermal screening length [ pO[ . 

To study the contributions from trapped and passing 
electrons to the charge density {—d'^cp/dx'^) and how such 
contributions are affected by various parameters, we show 
in Figure || plots of ntr and np and the charge density p as 
a function of x for several values of ip and S. Figure ^ (a) 
and (b) plot 100xnt,.(x), 100 x [np{x) — l], and 100xp{x) 
for {ip, S)~{2 X 10^^,0.1). For an ambient plasma with 
Te — 700ey and Ad = 100m as found at ionospheric 
heights by FAST satellite in the environment of BGK 
EHs |l^, this case corresponds to -0 = 1-4 x lO^^V and 
d — 10m. As shown, in this weakly nonlinear case, the 
maximum perturbation in Up is only 0.5% and in ntr 
0.4%. The perturbation in the charge density p is < 
0.2%, and occurs all within one Xjj. 

np{x), ntrix), and p(x) for (■0, 6) — {b, 4.4) in Figure 
§ (c) and (d). This choice corresponds to a point close 
to the curve ftr(w = —ip) = in Figure 2 and is an 
extremely nonlinear case. One can see that the total 
charge density perturbation goes ~ 10% negative and 
~ 25% positive, corresponding respectively to electron 
density enhancement and depletion. With similar for- 
mat, Figure ^ (e) and (f) plots a case with same 5 and 
■0 = 1 to illustrate the change in rip, ntr, and p of an EH 
with equal width but smaller amplitude. By locating this 
case in Figure |[ one notices that farther away from the 
ftr{w = —ip) = curve, the dip in ntr is filled up and the 
charge density perturbation only increases to 5% positive 
and 2% negative. 

These examples demonstrate how trapped electrons 
produce negative charge density perturbations and pass- 
ing electrons positive charge density perturbations owing 
to the decrease in their number density. It is always true 
that ntr > 0, since the number density cannot be nega- 
tive, and therefore trapped electrons always contribute to 
negative charge density regardless of the strength of the 
nonlinearity. This result disagrees with the picture that 
the positive core is due to a deficit of deeply trapped elec- 
trons, and that this positive core is screened out by the 
ambient electrons [|j . It is also different from the conclu- 
sion that the trapped electrons are screened out by the 
resonant or nonresonant passing electrons depending on 
the EH velocity in ||l|. 

We now turn to the issue of minimum size of EHs. 
From the illustrations of Figure |3[ one sees that the 
charge density variation is the net balance of the negative 
charge produced by trapped electrons and the positive 
charge density produced by a depletion of passing elec- 
trons inside the solitary potential. The trapped electrons 
must distribute and oscillate in such a way to yield the 
desired negative charge at the flanks to shield out the pos- 
itive core. The entire solitary object is a self-consistent 
and self-sustained object with zero total charge and does 
not require screening from the ambient plasma. Thus, 
the size of EHs are not restricted to be greater than A d ■ 
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In summary, wc demonstrated that the positive core of 
the EH is shielded exclusively by trapped electrons oscil- 
lating between their turning points and resulting in the 
excess negative charge. A BGK EH consists of an elec- 
tron density enhanced region and depletion region, and 
this means that the total charge for a BGK EH is zero. It 
thus follows that two separated EHs do not interact and 
the concept of negative mass is not needed. There does 
not exist a minimum size for BGK EHs since they do not 
rely on the thermal screening from the ambient plasma. 
These features are independent of the particular choice 
of potential profile and passing electron distribution, and 
also independent of the strength of nonlinearity, because 
the only principle upon which the arguments are built 
is conservation of charge (the continuity equation). The 
restriction to empty-centered EHs which was used previ- 
ously is relaxed to obtain a more general width-amplitude 
relation of an inequality form. 

Finally, note that the arguments and results we ob- 
tained for BGK electron solitary waves can be analo- 
gously applied to BGK ion solitary waves. 
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FIG. 1. Please see the text for explanations. 




I I I I I I 

1 2^3 4 5 

FIG. 2. the width-amplitude relation of BGK EHs that are 
not restricted to be empty-centered for a Gaussian potential 
and Maxwellian ambient electron distribution 
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FIG. 3. Trapped electron density (ntr), passing electron 
density (rip), and charge density (p) for {ip,S) = (2 x 10~^, 0.1) 
in (a) and (b), {tp,S) = (5,4.4) in (c) and (d), and 
{ip,6) = (1,4.4) in (e) and (f). 
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